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Abstract. - In synchrotron Mossbauer spectroscopy, the nuclear exciton polariton manifests itself 
in the lineshape of the spectra of nuclear forward scattering (NFS) Fourier-transformed from time 
domain to frequency domain. This lineshape is generally described by the convolution of two 
intensity factors. One of them is Lorentzian related to free decay. We derived the expressions 
for the second factor related to Frenkel exciton polariton effects at propagation of synchrotron 
radiation in Mossbauer media. Parameters of this Frenkelian shape depend on the spatial config- 
uration of Mossbauer media. In a layer of uniform thickness, this factor is found to be a simple 
hypergeometric function. Next, we consider the particles spread over a 2D surface or diluted in 
non-M6ssbauer media to exclude an overlap of ray shadows by different particles. Deconvolving 
the purely polaritonic component of linewidths is suggested as a simple procedure sharpening the 
experimental NFS spectra in frequency domain. The lineshapes in these sharpened spectra are 
theoretically expressed via the parameters of the particle size distributions (PSD). Then, these 
parameters are determined through least-squares fitting of the line shapes. 



1. Introduction 
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1.1. Synchrotron Mossbauer spectroscopy and 
lengths distribution. — Electronic and structural 
properties of materials were widely studied by Mossbauer 
spectroscopy in recent 50 years. Nuclear probe placed 
into deep structure scale interacts with its environment 
enabling to quantify the hyperfine interactions with 
unprecedented energy resolution. Cutting edge time 
domain Mossbauer spectroscopy using a synchrotron 
source has a potential of even twice better resolution. 
This is because the synchrotron Mossbauer spectra 
show no lincwidth contribution from radioactive source, 
and only the linewidth contribution of absorber sample 
remains. For this reason the spectra obtained in time 
domain and Fourier-transformed to frequency domain 
might have narrower linewidths compared to conventional 
energy spectra measured in velocity scale. 

Another attractive feature of the time-domain 
Mossbauer spectroscopy, to be explored in this work, 
consists in the possibility of studying the particulate 



matter. The point is that the length of the pathway 
through the resonant media is also reflected in nuclear 
resonant scattering pQ. The NFS theory [1] predicts in 
time domain the occurrence of the dynamic beats with 
periods depending on the lengths of radiation pathway. 
The pathway-dependent beats interfere with quantum 
beats originating from the quantum hyperfine splitting of 
nuclear levels. 

In this work, the possibility of obtaining the informa- 
tion on the particle size distribution (PSD) from NFS is ex- 
plored. Distribution of electronic and nuclear densities can 
be characterized by small-angle scattering (SAS) of x-rays 
and neutrons, respectively. In contrast, the unique fea- 
ture of the method proposed here consists in determining 
the chemically-specific densities of resonant nuclei. Ferro- 
or aniferromagnetic ordering sensed by Fe spins in mag- 
netic materials is probed using the 57 Fe nuclei. Therefore, 
our approach is relevant to various natural or synthetic 
magnetic compounds, doped by iron isotope, or contain- 
ing iron intrinsically. Particularly, in application to mul- 
tiphase systems, where at least one of the phases contains 
iron, there appears a unique chance of getting deeper in- 
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sight into the interfacial structure. Thus, the proposed 
method has the potential to make a substantial contribu- 
tion to experimental research both in multi- and single- 
phase systems. 

Consider the nuclide-embedding particles placed into 
the synchrotron beam in a way that excludes the ray shad- 
owing of one particle by another. The averaged over dis- 
tribution chord lengths determines the total intensity 
of resonant scattering in forward direction. 

In SAS, the distributions of nuclear or electronic densi- 
ties d n (C), d e (() are proportional to the second derivative 
of the SAS correlation function 70(7"): 



d n , e {r) 



ld 2 7oM 

c" 



dr 2 



(1) 



Here 7o(r) = r y{r)/j(0) and 7(f) is self-convolution of 
density distribution [2]. The proportionality factor l/( 
is the inverse mean length among all the particle chords 
distributed between and the maximum particle dimen- 
sion D: 



c = 



(d(C)d( 



(2) 



In our case, ( is the average thickness of the resonant 
media. The condition of "no shadowing" makes the chord 
lengths distribution (CLD) of randomly oriented parti- 
cles equivalent to the resonant-media thickness distribu- 
tion (RTD). In samples without mixing the resonant and 
non-resonant media, the RTD is quite analogous to the 
thickness distribution introduced previously in x-ray spec- 
troscopy for thickness-related corrections of x-ray absorp- 
tion spectra (XAS) [3]. While the thickness distribution 
determines the XAS and NFS spectra, the CLD is known 
to underlie, except SAS, a variety of properties of porous 
matter: conductivity, transport and relaxation. Investi- 
gating the CLD is subject matter of a special geometric 
sphere of knowledge [I] . 

1.2. Frequency-domain NFS spectra: the line- 
shapes. — In this work, we limit ourselves to the case of 
full coincidence between thickness and chord lengths dis- 
tributions (RTD=CLD). For example, this is the case of 
single layer of particles spread over a surface or Mossbauer 
particles well diluted in non-M6ssbauer media to avoid 
shadowing of one particle by another. This case includes 
not only the spherical particles of various particle size dis- 
tributions (PSD's), but also a network of wires or fila- 
ments viewed as discs in the plane perpendicular to the 
beam propagation direction. The relationships between 
the PSD, CLD and curves of SAS are known as solutions 
of the corresponding inverse scattering problems [2 . In 
a similar way, we will establish the relationship between 
PSD, CLD and the parameters of lineshape in the Fourier- 
transformed NFS spectra, such as shown in Fig. 1 (b). 

Our approach to the NFS data treatment in frequency 
domain is very different from the conventional one, that 
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Fig. 1: (a) Energy dependence of the NFS intensity in case 
of scatterer with the sextet of resonance lines of the transi- 
tion probabilities 3:2:1:1:2:3. Due to the multiple scattering 
each resonance shows the double-hump structure at large z. 
The spectrum was obtained by squaring the absolute value 
of the Fourier-transformed sum of Eq. (3). It was verified 
that exactly same spectra can be obtained using the formula 
(Eq. [7]) for squared transmission amplitude: |1 — R S (E)\ 2 = 

1 1 — exp ^— i E-Ei+iT g /2 ) I • H ere Lo is the linewidth ( 
0.097 mm/s in 57 Fe) 



i B __ Bi+i r /2 / 

'• The individual Mossbauer thicknesses 



[iiz for each nuclear transition are indicated. Total dimension- 
less thicknesses of 24, 48 and 96 are given by the sum over all 
transitions n tot z = E/ijz. (b): Frequency dependence of the 
NFS Fourier spectra, obtained using the same parameters of 
thickness and hyperfine interactions (magnetic hyperfine field 
i/hf=500 kOe, quadrupole splitting e = 0). 



has been known since the first observations of quantum 
beatings [5HZ]- m the frequency domain, we show that a 
single quantum beat can be disentangled from the complex 
sum of beatings at the condition of large enough hyperfine 
fields (i/hf ^ 200 kOe) . Fitting the lineshape for any indi- 
vidual line of the Fourier spectra permits to determine the 
RTD parameters even for complex distributions of thick- 
ness of the resonant media. 



3. Prerequisites 
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Fig. 2: NFS spectra of Fe2C>3 and GdFe03 fitted using the 
theoretical intensities derived from the Eqs. ((5| and ((6]). The 
goodness of fit is always best for a — 0. Except the experimen- 
tal spectra and fitted theoretical spectra for a = 0, the sim- 
ulated spectra for a = 1 nm//im (Fe2 0a) and a = 2 nm//im 
(GdFeOs) are also shown. The unenriched oxides Fe2C>3 and 
GdFeOa contain the density of the 57 Fe nuclei by 100 and 230 
times smaller than that in bcc metal of a— 57 Fe. 



Fig. 3: Energy spectra of the scattered radiation corresponding 
to approximate (a) and exact (b) solutions of the wave equation 
for radiation pulse propagating through resonant media. The 
squared amplitudes of radiation field vs. energy are calculated 
using the Eqs. ([5]) and © in (a) and (b), respectively. Note 
the coincidence of spectra in (a) and (b) for H^f > 200 kOe. 



3.1. The approximation of singular magnetic 
hyperflne fields. — Already the first studies of NFS 
|10j showed that a phenomenon of coupling between 
quantum beats and dynamical beats often takes place at 
large lengths of radiation pathways through the resonant 
media. This hybridization results from broad magnetic- 
hyperfine-field distributions that have a strong effect on 
the time evolution of the intensity of the forward scattered 
radiation. This is the case of magnetic alloys having 
essentially disordered atomic structure, such as Invar, or 
fully amorphous alloys and compounds. Distribution of 
quadrupole parameters in nanoparticles or in deformed 
crystals is another example of inhomogeneous broadening 
that gives rise to the coupling between quantum and 
dynamical beats. In what follows, we disregard these 
cases to focus on the case of singular value of hypcrfinc 
magnetic fields for all the 57 Fe nuclei involved. 



2.2. The approximation of large magnetic hy- 
perflne fields. — Even for a singular values of hyperfine 
fields the radiation scattered forward in samples of large 
thickness exhibits a broad energy distribution around each 
hyperfine transition, showing the so-called double-hump 
picture of a hyperfine transition (Fig.l, a). This picture 
arises from solving the wave equation for propagation of 
a radiation pulse through a resonant medium [10] . Closed 
analytic solutions of corresponding wave equations exist 
only for a single-line resonance. However, if we remain in 
the limit of large energy separation between the hyperfine 
transitions, there exists a good approximation to the ana- 
lytic solution even in the case of multiple resonances. This 
approximation is applicable when the energy separation 
between different transitions is large enough compared to 
natural linewidth [5MlO|. 

Then the approximate solution for the field amplitude 
of the radiation scattered in forward direction can be ex- 
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pressed as follows [TU] : 

E s {t,z) oc^^exp(-|A^-I) Jl( ^g? ) (3) 
j " * -\/Pi ZT 

Here the index Z numerates the nuclear transitions with 
the energy difference of Ai5; between ground and excited 
states. The factor e~ r / 2 describes the free-nucleus decay. 
The dimensionless time r = t/to is expressed in units of 
free decay time to (e.g., io=141.1 ns for 57 Fe). The nu- 
clear absorption coefficient \ii = cropfhMVWi is expressed 
individually for each transition using the weight factor wi 
of the Z-th transition. The factor wi takes into account 
the probability partition between nuclear sublevels split 
by hyperfine interactions, so that = 1. Other factors 
are the resonance cross section er , the density of the res- 
onant nuclei p, the Lamb-Mossbauer factor /lm 5 and the 
resonant nuclide isotope abundance 77. 

The time and space variables are entangled in the Eq. ([3]) 
via the complex argument of the Bessel function of first 
kind and order one: 

Using the variable = piz/A the power series of the 
Bessel argument ^fji\zf jl can be replaced by the power 
series of \f^\j ■ The dimensionless parameter individual 
for each Z— th transition £j is called the transition effective 
resonant thickness, resulting in radiation field: 

E*(t,z) oc J2^M-UE lt - T -) Jl{2 ^ ) (5) 

In terms of Eq.©, the spectra shown in Fig. 1 (a) 
represent the squared absolute values of the correspond- 
ing Fourier transforms \T{E s {t, z))\ 2 . The values of ef- 
fective resonant thickness \i\z are Z— specific and propor- 
tional to the relative intensity of the Mossbauer line for 
the Z-th transition. Three values of total Mossbauer thick- 
ness (fj, tot z = Yiipiz = 24,48 and 96) are used to draw 
the Fig.l. The magnetic hyperfine field of 500 kOe and 
zero quadrupole splitting are assumed. Very close hyper- 
fine parameters were observed previously in GdFeC>3 [IT] . 
Spectral weight of the quantum beat frequency spectra 
shown in Fig.l (b) represents the cos- Fourier transform 
of the squared sum in right-hand of Eq.©. Information 
contained in such Fourier images is the same as in experi- 
mental data. In what follows, we elaborate the concept of 
lineshape in these Fourier spectra. 

2.3. Are the shifts of quantum beat occuring in 
forward scattering?. — The (Eq© was introduced in 
1998 by Shvyd'ko et al (see, Eq. (16) in the work [TO]). 
In the Eq.©, the radiation field E s (t,z) exhibits at the 
time t = similar phase for each I th transition. The same 
assumption was used in Refs. [6], [HIIT3] . However, some 



of the previous works [7J[S] supposed the existence of a 
finite difference between the initial phases (at t = 0) for 
different transitions (see the Eq.(29) in Ref. [8]): 

l \ m / 

(6) 

Here tpi is the same factor as in Eq.©, tpi = 
Ji(2 v / 5r) / \f£jr and AEi m is the energy difference be- 
tween Z th and m th transitions. In the experimental data 
for magnetized 57 Fe foils, the finite quantum beat shifts 
were claimed in two early works [7J5] , but not reported in 
precedent [6] and succeeding [T0HI3] studies. Quite large 
values of the parameter a were reported for the NFS with a 
single QB component, a = 0.133 nm//im [7J and a = 0.14 
nm//xm [5] If such QB shifts were present for some of 
the QB components from sextet, the shape of the curve 
of sum-total NFS would distinctly depend on the thick- 
ness even in the absence of any distribution of thickness. 
In case of thickness distribution, the shift is expected to 
smear the quantum beats. On the other hand, no smear- 
ing is expected in case of validity of the zero-shift formula 
(Eq©, which accounts for the thickness distribution only 
through the envelopes rpi . 

No QB shifts were assumed previously at fitting the 
NFS data for powdered antiferromagnetic oxides [11] . Let 
us compare which of the Eqs. © and [6]) fits better the 
experimental time spectra. In Fig 2, the NFS spectra of 
Fe 2 03 and GdFe03 are fitted using both Eqs. © and 
[6]). The goodness of fit is always best for a = 0. Ex- 
cept the experimental spectra and theoretical spectra for 
a — 0, the spectra for a — 1 nm//im (Fe20a) and a = 2 
nm//im (GdFeOa) are also shown. Both these oxides were 
prepared with natural abundance of 57 Fe. The density 
of 57 Fe nuclei in such unenriched F82O3 and GdFe03 ox- 
ides is lower than in 57 Fe foil by the factors of 100 and 
230, respectively. For each pair of transitions, Z and m, 
the terms a£ m (r£ m / AEi m ) are relatively small. In case 
of magnetic sextet of transitions with small quadrupole 
shift, there are three groups of sums E m (r£ m / AEi m ) with 
relative values of ~ 0, ~ 1, and ~ 2. Each of these 
terms aE m (r^ m /A£'/ m ) is smaller than the QB shift for 
a = 0.133 nm/ /im in 57 Fe foil at least by an order of mag- 
nitude. Despite very small values of a assumed, Fig. 2 
shows strong divergence between experimental and simu- 
lated for a / spectra. Shape of the time dependence 
of the squared field amplitude E^(t,z) calculated from 
Eq.© exhibits very strong changes even for small devia- 
tions of a from 0. Such a variation of shape is the best 
measure for the value of a. No significant deviation of a 
from is observed in our data shown in Fig. 2. This is 
contrasting to the interpretations derived from the data 
in magnetized 57 Fe foil with just two unsuppressed tran- 
sitions [7j. The complex shape of the NFS time spectra 
for magnetic sextet is a more reliable observation than the 
shift of QB reported for the magnetized 57 Fe foil doublet. 
In fact, experimental NFS data suffer generally from the 
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uncertainty of time origin. This can be caused, for exam- 
ple, by a finite spread of prompt pulse. This uncertainty 
may lead to a misleading shift of spectra as a whole, how- 
ever, may not change the shape of the NFS pattern. 

Thickness distribution modeled by the exponential func- 
tion d(£) = ^q 1 exp(— CCo" 1 ) is consistent with the NFS 
spectra of Fe2C>3, but inconsistent with the NFS spectra 
of GdFeC>3 (Fig. 2). Except difference in quadrupole split- 
ting (cf. e = —0.1 mm/s in Fe2C>3 and e = +0.03 mm/s 
in GdFeC>3 [11]) these samples differ by factor 4 in their 
Mossbauer thickness and by factor of 10 in their metric 
thickness. The thicker sample GdFeO*3 showed a disagree- 
ment with the assumption of exponential distribution. In 
such a thick sample, the particles were most probably 
stacked in several layers, shadowing each other. Resulting 
distribution is much closer to (S-function than to exponen- 
tial distribution. We obtained £ ~ 22 assuming the uni- 
form thickness, = — £ ). By contrast, in the thin- 
ner sample ( Fe2 0*3), the distribution of thicknesses about 
the mean value £ was relatively broad. As the result, the 
fitting quality was equally good for d(£) = — £ ) and 
d(£) = ^q 1 exp(— £/£o) when the fitting is done in the time 
range between 30 and 130 ns (Fig .2). However, the ex- 
ponential distribution suits better when the spectra were 
fitted in broader range. The parameters can be fitted ei- 
ther in time domain or in Fourier-transformed spectra. 
Below we present fitting the Fourier spectra as the most 
convenient and transparent method. 

2.4. Exact energy spectra of the scattered ra- 
diation compared with the large-fields approxima- 
tion. — The condition of validity of the approximation 
(EqJSJ) was given in 1998 by Shvyd'ko et al [10]. It was 
pointed out [I0J that the Eq.© is valid even in the thick 
samples, if the nuclear transition energies Ei are well sep- 
arated from each other. This means that the condition of 
validity of Eq.(3) is the greatness of the magnetic hyperfine 
field compared to natural Lorentzian linewidth. Numer- 
ically, this statement is verified and confirmed in Fig. 3. 
Our calculations aimed to answer the question: how large 
fields arc sufficient to apply the approximation? In Fig. 3, 
the squared value of |1 — R S (E)\ is plotted, where R S (E) 
being the system response function: 

R s (t,z) =exp f -i V ^ -] (7) 

V ' ; V \ Z ^E-E l + ir /2 J V ; 

The Eq.JT]) is the exact result, therefore, the Fig. 3 
presents the comparison between the approximate result 
\T(E s (t,z))\ 2 (a) and the exact one |1 - R S (E)\ 2 (b). 
Large value of thickness fi tot z — 96 was chosen to plot 
the Fig. 3. Full coincidence of the exact result and ana- 
lytic solutions is shown for H^f > 200 kOc. 



3. Results 



3.1. Lineshapes of individual lines of the Fourier 
spectra for hyperfine sextet. — In large hyperfine 
fields (i?hf ^ 200 kOc), the difference between numeric 
and analytical solutions is negligible. Therefore, the 
squared sum over hyperfine transitions (Eq|3]) can be 
converted into sum over the intertransitional pairs: 

I(t, z) cx z 2 e~ T A n (z, t) cos zu n t (8) 

n 

Fourier transformation of the Eq.© results in the quan- 
tum beat frequency spectra shown in Fig.l (b). In case 
of pure magnetic hyperfine interactions, the quantum 
beats frequencies w n and amplitudes A n (z,r) can be 
expressed via the parameters O = jg3/2^NH\ 1 {/h and 
i>j = • / i( 2 \/jCminT)/\/jCminT (j = 1,2,3) , respectively, 
as shown in the Table 1. 



Table 1. The quantum beat frequencies vo n and am- 
plitudes A n {z,T) for nuclear levels split by purely mag- 
netic hyperfine interactions. The frequencies are expressed 
through O — \g3/2fJ'NHhf/h and the amplitudes are ex- 
pressed through ipj = J\ (y/jfizr) / \/JjIzT (J = 1,2,3). 
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There are eight spectral lines corresponding to eight 
nonidentical intertransitional energy differences for purely 
magnetic hyperfine interactions. According to the Table 
1, each line in the frequency domain has its own lineshape. 
The line at highest frequency is the widest one. This line 
turns out to be most suitable for our analysis because 
the Fourier transformation of 0| = J 2 / ^/^t = 

Jf (2-y/3£ m i n T)/3£ m i n r can be found analytically. Here we 
denoted by £ and fj,z the transition-specific quantities £ = 
3£ m in = /iz/4 = // ot z/16.By the minimum thickness, 
£min — M tot;z /48 we denoted the thickness of the weakest 
transition (1 in the notation 3:2:1:1:2:3). 

In presence of quadrupole interactions some of these 
eight lines split into doublets and triplets to produce four- 
teen lines, however, the line of highest frequency always 
stands alone. This outstanding feature makes it easy to 
apply the suggested here analysis in frequency domain. It 
is from this line the parameters of thickness distribution 
are derivable most easily. 
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Fig. 4: Thickness dependence of the full-width-at-half- 
maximum (FWHM) of the highest-frequency line shown 
in Fig. 1(b). The lineshape is the convolution of the 
thickness-independent Lorentzian and thickness-dependent 
term 1 1 - a F 3 [{§, |}, {§, 1, §}, with 2 F 3 being the hy- 

pergeometric function. Pure Lorentzian contribution to the 
linewidth at z = is 2.26 MHz. 



Fig. 5: Comparison of the universal distribution functions nor- 
malized by both area and height: the " Polaritonian" , ob- 
tained as Fourier transform of the QB envelope J-(ipi) = 
1 -3 F 2 [{\, f }, {§, 1, §}, - 4( £ 2 J o)2 ] (continuous line), com- 
pared to Lorentzian ^ i+^-Lq)' 1 /t' 1 (dashed line) and Gaus- 
sian . 1 exp(— ) (dotted line). 

V27TO" 2°" 



3.2. The shape and width of highest-frequency 
line. — In Eq. ((SJ, except the functions ipji the factor 
e~ r contributes the lineshapes. Unlike the ipj— dependent 
factors individual for each line, the factor e~ r is common 
for all NFS data, therefore, it can be easily removed by the 
sharpening procedure as shown below. Here we note that 
the e _T -factor dominates the full-width-at-half maximum 
for small linewidths (FWHM iS 5 MHz). Remaining fac- 
tors shown in the Table 1 become predominant for larger 
linewidths. The Fourier images of these remaining fac- 
tors are the A— shaped functions with some shoulders. In 
Fig. 1(b), the shoulders start to be visible for the singlet 
line as the thickness increases. In the limit of zero thick- 
ness the Lorentzian contribution to the linewidth is 2.2575 
MHz (Fig. 4). The plot of FWHM vs thickness is not lin- 
ear because the thickness-dependent component differs in 
shape from the Lorentzian L = J-{e~ T ). 

Since the NFS spectra are defined for t > the Fourier- 
transform is defined as J-(f) — / °° f(t) cos(ujt)dt. 
For the sake of normalization of the line intensities in 
Fourier spectra we introduce the new variable f that is 
another useful parameter related to Mossbauer thickness: 

7T f tot 7T 

C = n£ = —Liz = - — = — u tot z . (9) 
^ 4^ 4 16^ v ' 

The new parameter ( presents the measure of FWHM in 
the NFS Fourier spectra vs. angular frequency w. In con- 
trast, the parameter £ is the suitable measure of FWHM 
in the NFS Fourier spectra vs. linear frequency. Trans- 
forming the NFS time spectra, Eq. (|5J|, to the frequency 



domain one obtains a set of lines, whose lineshapes are 
represented by the convolution of the Fourier images of 
the free decay component and the thickness-dependent 
component. First component is the free-decay Lorentzian 
L = J r [exp(— t)] independent of the polaritonic param- 
eter £. Second is the polaritonic component related to 
propagation of radiation pulse through the resonant me- 
dia. Because of this convolution the full lineshape can 
be found only by numerical calculations. For a single- 
line spectrum this was already done [11]. The spectral 
lineshape was shown to evolve with increasing thickness 
from the Lorentzian profile to a different profile (so-called 
A-shaped [llj). However, no analytic expression was ob- 
tained yet for the thickness-dependent component of the 
lineshape. This solution will be found in the present work. 
The method of determination of the thickness distribution 
parameters simply follows from this analytic solution. 

In very thin samples, the same analysis could be ap- 
plicable to the line at w n = 3£7 shaped as simply as the 
Fourier-image of ipl, however, this line is just a tiny satel- 
lite of the strongest line of the NFS Fourier spectra (see 
Fig.l,b). That is why the operational range of thicknesses 
for this line is narrower than that for the line at w n = 19S1. 
In sharpened NFS Fourier spectra, we could manage to 
find the analytic solution for the spectral shape of these 
two lines, the narrowest one and the widest one, but not 
for the other lines. 

3.3 Shape of the polaritonic component for a 
uniform-thickness layer (foil). — Thus, since the 
sharpened lineshape is well-defined theoretically one can 
fit the sharpened data with a thickness-dependent analytic 
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Fie. 6: Center of the Polaritonian 1 — /13 , — ^-1 vs. (u — 

wo)/C (Eq- II 1 P in semilogarithmic scale. Maxima and minima 
in the oscillating behavior of the lineshape are indicated by up 
and down triangles, respectively. 



Fig. 7: The widths Auj and the heights A7 of the satellite 
shoulders in the polaritonic lineshapes shown in log-log scale. 
The numeration counter of shoulders TV goes from outwards to 
the center of the spectral line. Shown on the plot are the power 
exponents calculated using N = 13. 



expression. Prior the Fourier transformation the experi- 
mental spectra are multiplied by the factor exp(r). Then 
the NFS sum for a uniform thickness ( (EqU} will have 
the highest-frequency term as follows: 

/ Jr (2VC^A) \ 2 
I{t, C) exp(r) oc C V J cos 190t (10) 

y vCt/tt J 

This procedure allowed us to obtain in the frequency do- 
main the lineshape deconvoluted from the natural decay 
Lorentzian. Therefore, in the frequency domain, instead 
of the intensities I (to), calculated previously [IT], we op- 
erate now with the sharpened intensities I s h (lj, £) = 
J-[I(t, C) exp(r)]. The Fourier transform of the remain- 
ing lineshape factor in the right hand side of the Eq. (fT0|) 
has the functional form: 



T 



4C 2 



(ii) 



Here lj is the dimensionless angular frequency correspond- 
ing to the dimensionless time r. By h^(x) we denoted 

the hypergeometric function 2Fz[{\, §}, {|, 1, §}, x], us- 
ing the last but one argument of the function as subscript. 
The value £ = 1 in Eq. (fTT|) correspond to the normal- 
ized, both in area and height, distribution density (Fig. 
5). Similarly normalized Lorentzian and Gaussian distri- 
butions are [l + tt 2 (lj — ljq) 2 ~\ and exp [— tt(lj ~ ljo) 2 ] . 
The FWHM of the polaritonic function is slightly larger 
than FWHM of Lorentzian but smaller than that for Gaus- 
sian. To within this difference caused by the difference of 
the lineshape, the natural linewidth caused by the natural 



decay exp(— r) is equivalent to the thickness of Ccqv=r = ft 
or £eqv=r = L The origin of lj is located at the highest 
quantum beat frequency, hence, the positive and negative 
values of lj are measured relative to ljq = 190, that is not 
shown in Eq. fTTj) . The same notation for lj will be used 
below (lj "measuring" relative to ljq). 

The polaritonic lineshape has a very sharp tip in the 
center and the wings decaying similarly to the Lorentzian. 
It has the main shoulders I\ « 0.7io at lj\ « w$ ± 0.2£. A 
sequence of less-marked shoulders exists in closer proxim- 
ity to the resonance. These shoulders are formed by shal- 



low minima and maxima at 



and 



respectively, 



as shown in semilogarithmic scale in Fig. 6. The max- 



ima are located at 



lj ±2C/Ntt 2 , N = l,2..oo.The 



minima are given by the roots of the equation: 



2C/V 



(wo ± uj^ 1 



= (12) 



. x - w™ n and 
JS ax - IT n de- 



The width of the shoulders Alj^ = lj\ 
the heights of the satellite peaks Aijv — *n ~ x n 
crease with N approximately as the power laws iV -19 and 
TV -1 - 3 , respectively. Both curves Aljn vs. N and AIn 
vs. N show a small curvature in the log- log scale (Fig. 7). 
They merge each other both in position and in slope as 
N increases and the oscillations of I(lj) fade out with lj 
approaching ljq. 

The NFS radiation field and the nuclear currents are the 
polaritonic subsystems feeding each other as the excitation 
pulse propagates through the resonant media [THUS]. As 
the time elapses after the arrival of prompt pulse of syn- 
chrotron radiation the frequency of the energy exchange 
between these two subsystems at a fixed coordinate z de- 
creases. The nodes and antinodes change each other less 
and less frequently. The region of small frequencies in 
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the NFS Fourier spectra is dominated by the harmonics 
from large r of the NFS time spectra. Lorentzian shape in 
frequency domain correspond to exp(— r) in time domain 
and, vice versa, the Lorentzian's decay asymptotic t -2 
would produce around Wo the triangularly peaked shape 
exp(— \oj — ojq\). In Fig. 5, we observe that the peak is 
ever sharper than triangular. Taking the function ip^ (r) at 
the points of antinode maxima we observe that ?/>|(T a ) de- 
cays with r as t -3 / 2 , i.e., even slower than the Lorentzian's 
asymptotic r -2 . Indeed, one can check that the Fourier 
transformation of the function (1 + |r|^) _1 gives a very 
sharp peak, similar to that of Fig. 5, but free of the oscil- 
lations investigated in Figs 6 and 7. 

3.4 Deconvolving the polaritonic component us- 
ing the technique of sharpening of the experimen- 
tal Fourier spectra. — Sharpening the NFS Fourier 
spectra is shown in Fig. 8 for the data collected from 
unenriched Fe2C>3 with natural abundance of 57 Fe. The 
spectra were first fitted in the experimental time domain 
0.18 < r < 3 and then extrapolated to the short-delay 
region, r < 0.18. Missing original data for the t < 0.18 
caused by the prompt flash overshooting detector were 
thus restored. The expression for fit function was given 
previously (see Eq.(28) in Ref. dJ). Instead of 8 lines 
for e = there are actually 14 lines for the combined 
magnetic and quadrupole Hamiltonian. The fitted quan- 
tum beats parameters are Hhf = 517.5 kOe and e — —0.1 
mm/s. Due to the smallness of e the lines are grouped into 
2 symmetric triplets, 2 symmetric doublets and four sin- 
glets. The doublets and triplets are clearly distinguished 
in the sharpened spectra. 

From the viewpoint of theory the exp(— r)-sharpened 
spectrum in the middle panel of Fig. 8 must possess the 
purely polaritonic lineshape. The factor exp(r) permits us 
to isolate the polaritonic lineshape in its undistorted form 
predicted by the theory. In the middle panel of Fig. 8, 
the largest-frequency line at vj n — 1957 possess the largest 
polaritonic linewidth £ = 3£ min = 4.86. 

The sharpening technique is also very useful for the il- 
lustration of the relationship between the width of dif- 
ferent spectral lines. To the first approximation we can 
neglect the non-linearity in the plot of FWHM vs thick- 
ness in Fig. 4 and assume the additivity of the Lorentzian 
and polaritonic linewidths. Therefore, deliberately, the 
spectrum might be further sharpened using the exponent 
larger than exp(r) in Eq. (|10p . Such a heuristic sharpening 
is shown for the sharpening factor exp(3r/2) in the bottom 
panel of Fig. 8. The residual linewidth of the narrowest 
line at w n = 357 became very small (£ res — 0.05). Such an 
oversharpening is equivalent to subtraction of 37r/2 from 
the linewidth, while the theoretic exp(r)-sharpening cor- 
respond to subtraction of it. According to the Table 1, 
the polaritonic linewidths determined by the function -0| 
is three times larger than the polaritonic linewidths de- 
termined by the function ip^ . This relationship gives us 
the following equation for the polaritonic linewidths of the 
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Fig. 8: NFS spectra of the Fe2C>3 powder and their unsharp- 
ened and sharpened Fourrier transforms. Top panel data are 
from [IT]. The data are plotted vs. linear frequency, to be dis- 
tinguished from uj in Figs. 1-3. The datapoints in lower panel 
are obtained via fast Fourier transforms of NFS time spectra 
according to formula Ij = ±'E k I k e~ l '- 27rj/m ^ k with the number 
of channels in NFS spectra m=4096. Solid (blue) line is the 
sum of Lorentzians of different widths ^ H fitted individually for 
each of 14 lines. Each Lorentzian is the Fourier image of the 
time-domain exponential function fitting the individual quan- 

.The 



turn beat envelope: T ex p( — Coff 1 "/ 71 " 
constraints between the parameters of £ e ff for different lines 
was used to reduce the number of parameters from 14 to one 
thickness parameter (see Eq.28 in Ref. [ll]h 



lines at zu n = 357 and 1957: 



3(0.05 + |) 



4.86 



(13) 



This relationship is in full agreement with our observation 
in Fig. 8. 
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In the middle panel of Fig. 8, there appear only a few 
experimental points per FWHM of a spectral line. Two 
questions are remaining in this respect. First of all, it is 
natural to ask how the appearance of spectra can be im- 
proved. What kind of experimental limitations should be 
lifted over to raise the number of experimental points per 
FWHM? Second, we must indicate the way how to apply 
the analysis in practice for determination of the particle 
distribution characteristics. 

3.5 Effects of temporal range cut-off and reso- 
lution in frequency. — Generally, the interval Aw be- 
tween points of the discrete Fourier spectra is inverse total 
sample time: 



Aw 



2tt 



(14) 



In NFS, the value r max = toAt is limited by the interval 
between bunches of electrons of the storage ring. The data 
of Fig. 8 were collected in single-bunch regime [11] , with 
the interval between bunches of 500 ns and r max = 3.54. 
This experimental restriction on T max = mAr determines 
the number of points per FWHM. The number of chan- 
nels m = 4096 defines the NFS sampling frequency. Ac- 
cording to the 'uncertainty principle', Eq. ([T4|) . the in- 
crease in the value of m will not influence the number 
of points per FWHM. However, in very thick samples, ap- 
proximately for £ > 700, the existing number of channels 
m = 4096 would provide insufficient sampling frequency 
for the NFS signal. To obtain the full set of frequencies in 
the Fourier spectra, Nyquist-Shannon sampling theorem 
prescribes sampling the signal with the frequency twice 
larger than the highest frequency of the signal. In this case 
only, one should increase the number of channels compared 
to our value of 4096. 



Thus, in Fourier spectra, total sampling time r max lim- 
its the number of points per linewidth, which is propor- 
tional to £. With increasing £, the broader is a line of 
Fourier spectra, the narrower is the function (r) = 

^Ji 'Qr i 'it j I \J Ct /ir^j in time domain. In the exper- 
imental interval < t < r max , the number of nodes and 
antinodes of the ^|(r)-function increases with increasing 
£. Table 2 shows what should be the value of thick- 
ness £ to match the experimental window T max = 3.54 
with the n-th node or n-th antinode of the function 

2 



(ji (2VC^A) 



Table 2. Values of thickness C and exp(r)— weighted 

intensity n ^ J\ ^2-\/(" r / 7r ^ /C T f° r the NFS node and 
antinode locations at the upper boundary of the experi- 
mental time window t = 500 ns (r « 3.54) in a layer of 
uniform thickness. 



No. 


Node C 


Antinode £ 


Antinode 


t =0.5/xs 


t =0.5 /is 


I(t— 0.5/^s) 
/(t=0) 


1 


3.26 


5.85 


0.0175 


2 


10.9 


15.7 


4.16-10~ 3 


3 


23.0 


30.0 


1.60-10~ 3 


4 


39.4 


48.6 


0.78-10- 3 


5 


60.2 


71.6 


0.44-10- 3 



In Fig. 9, the function (j x h^/CrJ^) /VOT™) is 
Fourier-transformed to the frequency domain for the £ 
values from Table 2 after the data cutoff at r max = 3.54. 
Since these Fourier- images are plotted against (a; — ujq)/C 
they are all close to collapsing to the single master curve 
(cutoff- free), shown above in Fig. 4. 

Three useful types of behavior in the cutoff-affected data 
(Fig. 9) are of interest. First, the number of points per 
FWHM is inverse of ( as discussed above. Second, the dis- 
crepancy between the master curve and the cutoff- affected 
data culminates in the line center, and decreases towards 
the line wings showing an oscillating behavior. Third, as 
shown in the inset, the oscillation drops with increasing 
(lo — uq) more rapidly when the value of C from Table 2 
is at node than at antinode of the (r)-function. All 
these simulations were done assuming the homogeneous 
thickness. When the resonant media consist of particles 
the cutoff effects will be smaller. This is not counterin- 
tuitive because the thickness inhomogeneity would smear 
the oscillations of ?/;f (t). In the frequency domain, the in- 
homogeneity of thickness would produce the effects similar 
to cutoff, namely, the more smeared oscillations of V'K 7 ") 
produce the less sharp vertice in related spectral line. In 
what follows, this effect will be compared for several model 
thickness distributions. 



3.6 Fitting the thickness as linewidth in sharp- 
ened Fourier spectra. — Although our sample was 
made of particulate Fe2C>3, first we consider it forming a 
uniform layer. Within the small-thickness approximation 
each of the ipj function can be expanded into power series 
of thickness (see Eq. 28 in Ref. QT|), that simplifies fitting 
£ directly in time domain. The refinement of thickness pa- 
rameter in this way has resulted in £ tot = 12£ mm — 19.5 
(£ tot = 6.19). This total thickness is slightly larger than 
the value of £^ ot = 4.13 calculated from the sample weight 
of 55 mg/cm 2 . The difference results from the inhomo- 
geneity of the thickness. This inhomogeneity is taken into 
account below using two approximations: (i) large spher- 
ical micronic-size particles; (ii) smaller particles of Fe2C>3 
(nanoparticles) filling the space between the larger parti- 
cles of non-resonant media (MgO filler). 



3.7. Thickness distribution and lineshape for fil- 
aments or wires of the radius R. — In the particulate 
matter, owing to the distribution of the chord lengths the 
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Fig. 9: Lineshape profiles illustrating the influence of finiteness 
of the experimental time window (cutoff effects) on the line- 
shape in the NFS Fourier spectra. The full curve is the cutoff- 
free data similar that shown in Fig.3. Points correspond to 

the Fourier images of the function ^Ji ^2-\/C' r / 7I "j /•y/C T / 7r ) 
cutted off at r = 3.54 for the values of ( from 3 upper lines of 
Table 2. 



theoretical NFS lineshape acquires the forms shown in the 
last columns of Fig. 10. In first and second columns, 
the PSD and CLD are shown, respectively. The hyper- 
geometric function of Eq. (fTTj) is reproduced in (a) in the 
same raw with the uniform thickness layer PSD and CLD, 
both represented by (5-functions. The second raw (b) cor- 
respond to the filaments or wires of uniform radius. They 
are placed into the plane perpendicular to the beam di- 
rection, therefore, their CLD correspond to the distribu- 
tion of the chord lengths ( of the discs in 2D. For a disc 
of the size 2R using the running impact parameter b we 
can express ( = 2\/R 2 — 6 2 and b = \/R 2 — C 2 /4- Sim- 
ilarly to thickness distribution [3] the CLD d{Q can be 
found from the relationship d{Q — dS/S. Here dS is 
the relative part of the sample having thickness There- 
fore, the disc CLD is d{() = S^dS/df) = \db/dC, \ = 
C(AR)~ 1 (R 2 - C 2 /4)-2. The lineshape can be obtained 
through integration of the disc CLD with the right-hand 
side of the Eq. (II]) and ( from the Eq. (fit)]) : 



2R 



1-/13 - 



I ( n*u 3 ) 



C 2 d( 



(i? 2 -C 2 /4)5 4i? 



nR 



16i?' 



(15) 

The last but one argument of the hypergeometric func- 
tion is used as a subscript in our notation li2{x) = 
2F 3 [{if},fi,l,2},4 

3.8. Lineshape for spherical particles of the ra- 
dius R. — In the third row (c) and in rows below (d,e,f) 
of Fig. 10, the scattering matter consist of OD-objects, i.e., 




(ro-G) K 



Fig. 10: Size distributions (left column), thickness distributions 
(middle), and NFS lineshapes (right) for 2D layer (a), ID wires 
(b), uniformly-sized spheres (c), spheres with rectangular (d), 
triangular (e), exponential (f) size distributions and with the 
size distribution of the form F(r) = (R/r 2 + 2/r) exp(—2r/R) 
(g). In (a) R is the layer half-thickness, elsewhere R is the 
average radius of sphere or wire. 



spherical particles. In (c) all the particles have spheri- 
cal shape and constant radius. The RTD function coin- 
cides with g?(£) when the particles arc distributed over 
the plane surface S in such a way that shadowing of one 
particle by another is excluded. Again, we have the im- 
pact distance b = \J R 2 — C 2 /4, however, dS/S is now 
{-KR 2 )- 1 {2Mb/dC) = C/2i? 2 (0 < C < 2i?). Then 
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4C 2 



C 2 d( 
2R 2 



AR 
IT 



1 - ft 6 



16i?' 



(16) 

Similarly to Eqs. (|llll5[) the last but one argument is used 
as a subscript hs(x) = 2 F 3[{|, |}, {|, 1, §},#]• 

3.9. Lineshapes generated by several charac- 
teristic particle size distributions. Assuming the 
spherical particles to have some distribution in size we 
observe (Fig. 10 (d, c) that the PSD's of rectangular 
and triangular shape result in the lineshapes with smeared 
shoulders. No shoulders can be apperceived for the cases 
(f) and (g), related to exponential PSD and exponential 
CLD, respectively. All the PSD's except (g) are normal- 
ized to have the same area and the same mean particle 
size (r) = R. The CLD's d(() are derived from PSD's 
F(r) using the formula [IB] : 



oo 

m = / F <r) dr 



(17) 



C/2 



In case of exponential PSD (f), i.e., F(r) — 
(l/R)exp(-r/R), the Eq.flJU) results in the CLD d(() = 
(C/4i? 2 ) exp(-C/2i?). Integrating this d(C) with the right- 
hand side of the Eq. (fTTj) and ( from the Eq. (|T0l) we 
obtain the spectral lineshape 



/ 



AC 



C 2 exp(-C/2fi)dC 
AR 2 
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20Fi?sin(3$/2) ^cos(<&/2) 
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64R 2 \i 



64 J? 2 



(18) 
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Fig. 11: The section from 170 MHz to 200 MHz of the quantum 
beat frequency spectra (square-shaped symbols) fitted with the 
right-hand expressions in Eqs. (|18p and (|19p (lines). 



The absorbers for NFS are frequently prepared through 
mixing the resonant nanoparticles (e.g., 57 Fe203 in |11| ) 
with larger particles of non-resonant filling medium. 

The divergence of F(r) at small particle radius (r — ► 0) 
is so strong that the integral J F(r)dr is also divergent. 
This nonitegrability implies that the spherical shape as- 
sumption (Eq. I17[) is rather unphysical for the exponen- 
tial CLD d(C) = (l/.R)exp(-C/-R) [IS]. Indeed, this d{Q 
is the chord length distribution in the background mate- 
rial filling the interstices between spheres [TB] ■ Integrating 
this d(() with the right-hand side of the Eq. (fTTj) and ( 
from the Eq. (fTUl) we obtain the lineshape of such a back- 
ground: 



Here $ = arctan(8i?/7ru;). As far as the relative weight 
of smaller particles increases from top to bottom of Fig. 
10 the shoulders became less articulated and disappear in 
(f)- 

Further increase in the fraction of small particles is the 
case (g) related to divergence of F(r) at r — > 0. Substi- 
tuting F(r) = {r/R 2 +2/r) exp(-2r/i?) into Eq. [J7]we ob- 
tain the simple exponential CLD d(() — R^ 1 exp(— £/2i?), 
which is broadly found among the natural particulate ma- 
terials. Since this F(r) has a singularity at r — > , 
such a PSD cannot be normalized in the same way as the 
PSD's (d), (e) and (f). However, d(() can be normal- 
ized in a usual way J d(r)dr = 1 with the average chord 
length (0 — R- Many kinds of complex interfacial me- 
dia exhibit the exponential CLD. Numerous examples of 
the exponential CLD's were demonstrated [T7] among the 
porous borosilicate glasses, heterogeneous catalysts, ce- 
ments etc. Also, if we mix the hard spherical particles of 
non-resonant media with soft resonant background we ex- 
pect an exponential CLD for the resonant component |18j . 



1 - hi 



A( 2 



e~iz(d( 



R 



R 



/ttcos((/>/2) 



(tt 2 



lli/,'- 



Here tj> — arctan(4i?/7rw). 



(19) 



3.10. Determination of the particle size distri- 
bution parameters. - The Eqs. (JTTJ), (HSJ), (O and 
(|19[) were employed to fit the upper-frequency line from 
the exp(r)-sharpened spectrum of Fig. 8. First, assuming 
a layer of uniform thickness the Eq. ([TT| was used. Fit- 
ting the line at 185 MHz with the hypergeometric function 
1 — h 3 / 2 (—A( 2 /it 2 oj 2 ^ has resulted in the value of the di- 
mensionless thickness Co = 4.36 (Fig. 11). According to the 
Eq.© the total Mossbauer thickness is C($ ot = 4(o - 17.4. 
This is very close to the value of C tot = 19.5 obtained 
above in time domain through fitting the NFS spectrum 
as a whole. 

Turning from dimcnsionless £ tot to metric units of z = 
C tot / 7r M tot and using /i tot = ctoP/lm?? with a — 2.56 • 
3.96 • 10 22 cm" 3 , / LM = 0.8, 77 = 0.0212, 
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Fig. 12: Chord length distributions d(z) used to fit the highest- 
frequency line of the NFS Fourier spectra of Fig. 8 at the values 
of the dimensionless parameters: £o = 4.36 in d(£) = <5(£ — £o), 



R = 2.89 in d(£) = 
R = 0.82 in d(C) 



c 

t C 



R = 2.38 in d(() = jre'R and 



we obtain the total coefficient of nuclear absorption for 
unenriched hematite /i tot = 0.172 /im _1 . From the Eq.Q 
the layer thickness zq that corresponds to Cq ' = 17.4 is 
(shown by the position of (^-function in Fig. 12): 



4C<$°77r A i t ° t 



16Co/7iV ot = 129/im (20) 



Certainly, a sample formed by the particles of Fe2C>3 
of a finite size appears worthy of a better characteriza- 
tion. Assuming all the particles having the same size and 
spheric shape, we can estimate the radius of particle us- 
ing Ea. ([15]). The fit results in the value of R of 2.84. In 
this model, the thickness is distributed linearly between 
and 2R (Fig. 12). The particle diameter of 168 /urn is 
larger by 30% than uniform thickness zo from the Eq. (f20|) . 
Here again as in Eq. ([20l) the hematite-specific length fac- 
tor r — 16 /nfi tot — 29.61/im brings the diameter 2Rr to 
the metric scale. The goodness of fit becomes better when 
we proceed from uniform layer to the ball-shaped parti- 
cles, namely, the coefficient of determination R 2 increases 
from 0.991 to 0.994. Although the increase is not large, 
the usage of the R 2 — coefficient is justified at comparing 
the goodness of fit for different models because all of our 
models use the same number of parameters. 

Among the models with distributions of particle sizes 
two models were tested, according to the Eqs. (|18|19|) . 
Goodness of fit was slightly better for the first (Eq. [T51 
R 2 = 0.995) than for the second (Eq. [TH R 2 = 0.993). 
Corresponding chord lengths distributions are also shown 
in Fig. 12. The quality of fit between these two distribu- 
tions is not much different, the values of fitted radius R 
for the CLD of spheres is nearly 3 times smaller than the 
value of R for the CLD of the background between spheres. 
In the chord length range 38 /im < z < 260 fim the CLD 
density for spheres is larger than the CLD density for the 



background. Opposite is true elsewhere. Clearly, the ex- 
cess in the range 38 /im < z < 260 /im is compensated 
by the contribution of large particles, although their CLD 
density above 260 /im is below 5%. In this sense, the chord 
length of 260 /im is the characteristic length invariant of 
the fitting model employed. Interestingly, this "inhomo- 
geneous" length is twice larger than the value obtained in 
Eq. (|20"|) for the homogeneous layer. 



4. Conclusions. — Starting from several charac- 
teristic states of scattering matter (foil, filaments, 
particles) we have solved the problem of determination 
of corresponding spectral lineshape parameters for the 
nuclear forward scattering Fourier spectra. The method 
is useful in the range of thicknesses where quantum beat 
envelopes vary slowly compared to the quantum beats 
themselves. In the frequency domain, we are able to 
predict the lineshapes determined by Fourier transforms 
of the envelopes. These lineshapes are of polaritonic 
nature since they originate from the energy exchange 
between the radiation field and nuclear excitation. Both 
the radiation and the nuclear currents are the compo- 
nents of the compound quasiparticle termed the nuclear 
exciton polariton. From these lineshapes we derive the 
parameters of resonant thickness distribution. 

When fragments of resonant and nonresonant media are 
interleaved so that the standard methods of neutron or 
x-ray small-angle scattering (SAS) cannot distinguish be- 
tween fragments of similar electronic or nuclear density, 
the proposed technique would be able to resolve between 
them. Therefore, if a fraction of the 57 Fe-containing 'par- 
ticles' cannot be detached from a 'membrane', then the 
methods of SAS would produce the information on the 
density distribution in the system 'particles+membrane'. 
Separately, the particles subsystem can be studied by NFS 
without membrane detaching that is most frequently un- 
realizable. Supported catalysts present the archetype of 
future applications. 

Practical applications of the proposed analysis are quite 
feasible already at the existing beamlines of the syn- 
chrotron rings of third generation. The sampling fre- 
quency does not require any improvements, because it 
can be needed only for very thick samples. The time res- 
olution of the avalanche photodiode detectors (~0.1 ns) 
correspond closely to the currently applied sampling fre- 
quency. Much more important for the proposed technique 
would be the improvements in the high-rate characteris- 
tics of the detectors and in the span of the time window 
between electron bunches. Change from multiple-bunch 
to single-bunch regime is crucial to achieve a better cutoff 
time T max . Also, with the advent of novel avalanche de- 
tectors and detector arrays [19l[20] the region of missing 
data at small r might be shortened. 

Inverse problem of finding the characteristic parameters 
of the particle chord length distribution from the polari- 
tonic lineshapes would be possible to solve using regular- 



-12 



Nuclear forward scattering in particulate matter: dependence of lineshape on particle size distribution 



ization methods. Then, not only PSD parameters could 
be refined starting from one or another hypothesis, but 
also the full PSD profiles could be reconstructed. 

This work was supported by RFBR-JSPS Grant 07-02- 
91201. 
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